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Abstract 

To resolve various outstanding issues associated with the twist four longitu¬ 
dinal structure function F£ =i (x) we perform an analysis based on the BJL ex¬ 
pansion for the forward virtual photon-hadron Compton scattering amplitude 
and equal (light-front) time current algebra. Using the Fock space expansion 
for states and operators, we evaluate the twist four longitudinal structure 
function for dressed quark and gluon targets in perturbation theory. With 
the help of a new sum rule which we have derived recently we show that the 
quadratic and logarithmic divergences generated in the bare theory are related 
to corresponding mass shifts in old-fashioned light-front perturbation theory. 

We present numerical results for the Fi and Fj j structure functions for the 
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meson in two-dimensional QCD in the one pair approximation. We discuss 
the relevance of our results for the problem of the partitioning of hadron mass 
in QCD. 


I. INTRODUCTION 

An important problem in applying QCD to deep inelastic scattering is the existence 
of power corrections to scaling, more commonly known as higher twist effects. They are 
essential for making precision tests of QCD. From the early days of the establishment of QCD 
as the underlying theory of strong interactions, the importance of a proper understanding of 
power corrections was recognized [|l] . Subsequently, leading ^2 corrections to the unpolarized 
leading twist structure function F 2 and the longitudinal structure function F L were analyzed 
by the Operator Product Expansion (OPE) (2J and Feynman diagram [pj] approaches. Later 
Qiu H has given an alternate method based on special propagators (utilizing unique features 
of light-cone coordinates) to simplify the analysis. 

Power correction to F L is especially interesting since the leading (twist two) contribution 
to F l is perturbative in origin in contrast to the case of F 2 . Thus the first, non-perturbative 
contributions to Fl occur at ^2 order. The complexity of the problem of higher twist appears 
in the OPE analysis which utilizes a collinear basis, since at twist four, there appears a 
proliferation of operator structures. In the Feynman diagram approach it has been shown, 
using a transverse basis, that contact could be made with light-front current algebra analysis 
with the result that the twist four part of Fl is given by the Fourier transform of the hadron 
matrix element of the minus component of the bilocal vector current. Since the minus 
component of the current involves the constrained fermion held, the relevant operator has 
explicit dependence on the interaction in contrast to the well-known result for the leading 
twist contribution to F 2 which involves the plus component of the bilocal current. Even 
after many years of investigation, an intuitive physical understanding of the interaction 
dependence in the structure of Fl has been elusive. We provide herein a physically intuitive 
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picture. 

Another important problem of current interest is the perturbative aspects of the twist 
four matrix element. Simple power counting indicates that in the bare theory the twist four 
matrix element will be afflicted with quadratic divergences [fj. Understanding the origin 
and the nature of these divergences will be quite helpful in finding procedures to remove 
them ( the process of renormalization). 

A third motivation to study the twist four part of Fl comes from the present status of 
deep inelastic scattering experiments. Measurements § of the ratio of longitudinal to trans¬ 
verse cross section in unpolarized deep inelastic scattering show 0 that power corrections 
play an important role in nucleon structure experiments in the SLAC Kinematic range. It 
is important to go beyond phenomenological parameterizations for a proper understanding 
of the non-perturbative nature of these corrections. 

Light-front analysis of deep inelastic scattering provides an intuitive physical picture of 
various structure functions at the twist two level. Recently, the resolution of an ambiguity 
at the operator level and the parton interpretation of the transverse component of the 
bilocal current have been achieved in an approach based on light-front held theory [p]. The 
physical picture of the transverse polarized structure function |J and a critical examination 
of the Wandura-Wilczek sum rule in perturbation theory [jllj have also been provided in 
the same approach. Both nonperturbative and perturbative issues can be addressed in the 
same language in this formalism which uses the Fock space expansion for all the operators 


and multiparton wave functions for the state |lT|. The approach also provides insights into 
various renormalization issues associated with the different components of currents and the 
Hamiltonian. 

In this work, we show that using the same framework, one can resolve outstanding 
issues associated with the twist four contributions to the longitudinal structure function. A 
brief summary of some of our results is presented in Ref. [[1^]. In this work we extend our 
previous calculations and also present several new results. Our starting point is the Bjorken- 
Johnson-Low (BJL) expansion for the forward virtual photon-hadron Compton scattering 
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amplitude. This leads us to the commutator of currents which we present in detail for 
arbitrary flavors in SU(3). Next we consider the specific case of electromagnetic currents 
and arrive at expressions for the twist two part of F 2 and twist four part of Fl in terms of 
specific flavor dependent form factors. In the rest of the paper we consider the flavor singlet 
part of the structure functions. We identify the integral of Fl ^ with the fermionic part of 
the light-front QCD Hamiltonian density. The consideration of mixing in the flavor singlet 
channel leads us to the the definition of the twist four longitudinal gluon structure function 
and then we find a sum rule, free from radiative corrections. 

The sum rule which the physical structure function has to satisfy involves the physical 
mass of the hadron which is a finite quantity. A theoretical evaluation of the sum rule which 
starts with the bare theory, on the other hand, will be afflicted with various divergences 
(see Sec. IV.) depending on the regulator employed. In order to compare with the phys¬ 
ical answer resulting from the measurement, we need to renormalize the result by adding 
counterterms. For the dressed parton target, for example, these counterterms are dictated 
by mass counterterms in the light-front Hamiltonian perturbation theory. For a dressed 
gluon target, calculations in Sec. IV B show that quadratic divergences are generated and 
one does not automatically get the result expected for a massless target. The divergence 
generated is shown to be directly related to the gluon mass shift in old fashioned pertur¬ 
bation theory. To a given order in perturbation theory, counterterms have to be added to 
the calculated structure function. The precise selection of counterterms is dictated entirely 
by the regularization and renormalization of the light-front QCD Hamiltonian. The choice 
of counterterms in the Hamiltonian, in turn, determines the counterterms to be added to 
the longitudinal structure function which results in a theoretical prediction of the physical 
longitudinal structure function. Recall that in Hamiltonian perturbation theory we cannot 
automatically generate a massless gluon by clever choice of regulators. The point we em¬ 
phasize is that the twist four longitudinal structure function is one to one related to the 
Hamiltonian density and that there is no arbitrary freedom in this relationship. 

We also note that in the pre-QCD era, there have been discussions about a possible S(x) 
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function contribution to the longitudinal structure function which may appear to invalidate 
the sum rule derived ignoring such subtleties. In two-dimensional QCD Burkardt has shown 
I3| that has a delta function contribution and he has discussed implications of this for 


the sum rule for pf. Obviously, will not be affected by such a singular contribution and 
we show explicitly in Sec. V that the sum rule is verified in two-dimensional QCD by virtue 
of the ’t Hooft equation. 

To gain an understanding of the nature of quadratic divergences, we evaluate the twist 
four longitudinal structure functions for quark and gluon target each dressed through lowest 
order in perturbation theory. The sum rule allows us to relate these divergences to quark 
and gluon mass corrections in QCD in time-ordered light-front perturbation theory. We also 
verify the sum rule in a non-perturbative context in two-dimensional QCD. We also present 
numerical results for F 2 and F L structure functions in this model using wave functions 
calculated in a variational approximation. Finally we discuss the relevance of our results for 
the problem of the partitioning of hadron masses in QCD. 

The plan of this paper is as follows. In Sec. II we derive the expressions for the twist 
two structure function F 2 and the twist four longitudinal structure function Fl using the 
BJL expansion and equal time (x + ) current algebra. The sum rule for F L is given in Sec. 
III. In Sec. IV we evaluate Fl for quark and gluon targets dressed through lowest order in 
perturbation theory and explicitly verify the sum rule. The sum rule is verified explicitly in 
a non-perturbative context in two dimensional QCD in Sec. V. In this section, to provide a 
qualitative picture, we also present numerical results for the F 2 and Fl structure functions 
in this model. In Sec. VI we discuss the issue of the breakup of hadron mass in QCD in 
the context of our sum rule. Discussion and conclusions are presented in Sec. VII and our 
notations and conventions are summarized in an appendix. 
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II. PRELIMINARIES 


In this section we present the expressions for structure functions for arbitrary flavors 
in SU(3) which follow from the use of the Bjorken-Johnson-Low expansion and light-front 
current algebra. In terms of the flavor current J^(x) = the hadron tensor 

relevant for deep inelastic scattering is given by 

wg = ^(PIWK), Jf(0)]|P). (2.1) 

The forward virtual photon-hadron Compton scattering amplitude is given by 

IT = */dV»yp|r(j£K)jr(0))|P). (2.2) 


We have 


T^b ( x ) Q 2 ) = 2 /°° dq' + 

J —OO 


q' + - q+ 


Using the BJL expansion flUfl , we have, 

1 


ryn/i V _ 

^ ab 


dcd^qp p*©, juo)] c+ .„ip) + 


(2.3) 


(2.4) 


where ... represents higher order terms in the expansion which we ignore in the following. 
In the limit of large q ~, from Eq. (Al), we have 




tjH— _ ^ p I (p±\2^{ab) , P -Q 

w ab - 2 F ^ b ) + 


* 2(ab )j 


(2.5) 


with x = -ff- and z/ = P.q. On the other hand, from Eq. 


Limit,-^ T+- = -- /dr^UP’^Pliy (0, 4"(0)]«+=o|P) ■ ( 2 . 6 ) 

q J 

The components of the flavor current Jf(x ) obey the equal- x + canonical commutation 
relation (to be specific, we consider SU(3) of flavors) 


Ja( x ), J b ( V) 


A 


= 2 ifabc Il>(x) 7 ^(x) S 2 (x ± — y^) S(x - y ) 


=y+ 


1 <9+ 

2 1 


e{x - y ) if abc V c (x \ y) + id abc V c {x \ y) 8 2 (x 1 - - y^) 


+ 2 ifabc e ( x “ y ) 


+ - id abc e(x - y ) d l x 


S 2 (x ± -y ± ) V c (x\y) - e” A J c (x \ y) 


V A J t 


S 2 (x ± - y L ) V c {x I y) + ^(z | y) 


(2.7) 
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In deriving the above relations, use has been made of the relation 


ifabc^c d a 5 c A c . 


( 2 , 8 ) 


We have defined the bilocal currents as follows. 


Vc(x I y) 

K(x I y) 
A£{x | y) 
~K{x | y) 


1 X-r, . A, 




2 + ^{y)— ^(x) 

-^(2/)y7^0) 

\ [^(aOyTWfa) + Tpiy) y7 /i 7 5 '0(^) 
^[^(a0y7 / VV(2/) ~ ^(y)^^ 5 ^) 


Further, we introduce the bilocal form factors 


(2,9) 


{/> i v"K i o) i p>= p“v^e, P4)+ev?(t 2 , p-o 
{P I VcK I o) I P)= P“vl(e, P4) + evl(4 2 , P4) 


( 2 . 10 ) 

( 2 . 11 ) 


^From Eqs. ([2.6|) and ( |2.7|) , we get, 


Limitq --too q T a b ‘2if abc P r c 


+ 2 J e ^A fabc < P I I °) I P ) + < P I V c (? I 0) I P) 

~ j «-ei« + f-e(r)[/^ <f I KK | 0) P) + cU (P | V c ({ | 0) | P) 


( 2 . 12 ) 


Note that matrix elements of A%{x \ y) do not contribute to unpolarized scattering. Using 
the dispersion relation given in Eq. Q2.3Q , together with Eqs. ( |2.5|) and Q2.12 ) and comparing 
the coefficient of q l on both sides, we get 


= -L-Jdi]e "’’[UcV'M + <t*V*(i;) 


(2.13) 


Comparing the coefficients of q + on both sides, we get, 


F 


L(ab) 


[X) = 


1 i ( q + ) 2 
P+ 

{P _ 

Q 2 7I-P+ 


>_L\2 


f d V e^ x [f abc (P | V“(£ | 0) | P) + d abc (P | V c (e I 0) | P) 

■ 2 J drje-^ x [f abc (P | V+(£ | 0) | P)+d abc (P | V c + (£ | 0) | P)]. (2.14) 
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We have introduced rj = bP + £, . 


Note that our result for Fl differs from the one given in the literature |1{J. The difference 
can be traced to the expression for Fl that one employs. It is customary [15|.|3[] to ignore 
target mass M 2 in the expression for Fl (see Eq. (|772|)). This leads to an incorrect expression 
for Fl which in turn will lead to an incorrect sum rule (see the following section). 

The electromagnetic current 


•I'M = JS(x) + yf-'s'M- 


(2.15) 


^Frorri the flavor structure of electromagnetic current, we observe that, only d a bc contributes 
to the structure functions in deep-inelastic electron-hadron scattering. Explicitly, we have, 

F 2 (x ) i 


x 


2nP + 


dr]e~ ir}X (P | V (f | 0) | P). 


(2.16) 


The longitudinal structure function is given by 
2 i (q + ) 2 6 


FAx)= — - 


Q 2 tt P+ 


d7]e~ ulx (P | V (e | 0) | P) 


i „2 


Q 2 7 lP+ 


x 2 / drje~ l71x {P I V (£ | 0) | P). 


(2.17) 


We have defined the functions 


10) = | 0) + ivtK I 0) + I 0). 


(2.18) 


In arriving at our final results we have used explicit values of the structure constants of 
SU( 3), 


11 2 
^338 — -^=, d$88 = ~^=, ^330 = d 880 = y g- 


(2.19) 


Vo is the flavor singlet component of the fermion bilocal vector current. 


III. SUM RULE 


Consider the flavor singlet part of the structure functions F 2 (f) and Fl(j) defined by 


8 
















F- 


2 (/) 


>X 


1 


X 


47 tP+ 


dye iVX (P | ['0(O7 + V’(O) -'0(O)7 + '0(Ol | P). 


(3.1) 


F f s 11 (g + ) 2 
L(/) ^ Q 2 n P+ 


drje l1)X {P | [-0(07 0(0) - 0(0)7 0(f) J I P) 

(P L ) 2 1 


: x 2 dije lvx (P | 0(O7 + 0(O) - 0(O)7 + 0(O | P). (3.2) 


Q 2 nP+ 

^From Eq. (ED it follows that F 2 ^f){—x) = F 2 (f)(x) and from Eq. (|3.2| ) we explicitly 
find that FJ^{—x) = —F[^(x). It can be verified |L2| that F[^ satisfies the sum rule, 


dx 


F T LU)M 2 ) 


x 


2 r (P | 9+-(0) I P) - FF(P I 7+(0) I P) 


Q 21 


(P 


(3.3) 


where 6+~ = i07 _ d + 0 is the light-front QCD Hamiltonian density and 6+ + = f07 + c) + 0 is 
the light-front longitudinal momentum density in the light-front gauge A + = 0. 

Here we have used the fact that the physical structure function vanishes for x > 1. 
Neglect of M 2 in the expression (Eq. (|A2|) ) for F L will lead to (P 4 -) 2 + M 2 instead of (P 4 -) 2 
in the above equation which would spoil the correct sum rule given below. 

pT =4 

The integral of L £ f) is therefore related to the hadron matrix element of the (gauge 
invariant) fermionic part of the light-front Hamiltonian density. This result manifests the 
physical content and the non-perturbative nature of the twist-four part of the longitudinal 
structure function. 

The fermionic operator matrix elements appearing in Eq. (|3.3|) change with Q 2 as a result 
of the mixing of quark and gluon operators in QCD under renormalization. Analyzing the 


operator mixing we obtain a new sum rule at the twist four level |12 

M 2 

-t'\- = 4 

/ o x 


rl ‘±Fr‘ = 4— 
1 Q 2 


(3.4) 


Where M is the invariant mass of the hadron and FI 4 = F[^ + P0 ff 4 , F ^ g j is the twist 
four longitudinal gluon structure function which we define as, 




1 xP + 


- e ~ : hP + y x 


dy e 


(P | (-)F +Xa (y~)F- Xa ( 0) + -g + ~F Xaa (y-)F\ tTa {Q) \ P) + (y 


0) 


(P 


_L\2 


(P 


472 


(P | (-)F +Xa (y-)F + Xa ( 0) | P) + ( y~ <- 0) 


(3.5) 
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WhereF MAa = d^A ua — d u + gf abc A^A v c . Note that in the definition of F£^(x) the second 
term where the arguments of F XrTa are interchanged is missing in Ref. [ |T2| . 

To our knowledge, this is the first sum rule at the twist four level for deep inelastic 
scattering or for QCD in general. The previously known sum rules in deep inelastic scattering 
are all at the twist two level. The operators involved are kinematical (light-front longitudinal 
momentum, light-front hclicity, etc.) in nature. In contrast, the sum rule we have derived 
involves a dynamical operator (light-front QCD Hamiltonian) thus revealing a new aspect 
of the underlying non-perturbative dynamics. Our results show that the measuremnent of 
the flavor singlet part of the fermionic contributions to the twist four longitudinal structure 
function in deep inelastic scattering directly reveals the hadron expectation value of the 
fermionic part of the light-front QCD Hamiltonian density in light-front gauge. 


IV. DRESSED PARTON CALCULATIONS 


A. Dressed quark with non-zero mass 


Next, we investigate the implications of Eq. (|3.3|) for quadratic divergences in 


m 


perturbation theory. We select the target to be a dressed quark and evaluate the structure 
functions to order g 2 . That is, we take the state | P) to be a dressed quark consisting of 
bare states of a quark and a quark plus a gluon: 


P,a)= b\P,a) | 0) 

dkfd 2 k 0 


£ 

ol,A2 


dkidHi ^f(P - k, - k,) 


2(2t r) 3 kt J y/2(27r) 3 AiJ 

02(-P, o | ki,a 1 -,k 2 ,X2)b\k 1 ,a l )a\k 2 ,\2) | 0). 


(4.1) 


In the previous work 0 we have given results for massless quark state. We have shown 
that the twist four longitudinal structure function has quadratic divergences in perturba¬ 
tion theory. In this section, we show that for a massive quark, in addition to quadratic 
divergences, logarithmic divergences are generated. We have, 
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Fl — Mi + M2 


(4.2) 


where 


Mi— 


dy-e~^ p+v X {P | \y~) 


ttQ 2 




^.[id 1 - + gA A 


+ 7 °m ^.[id 1 - + pA ± (0)] + 7 °m -0 + (O) + h.c. \ P), (4.3) 


and 


M 2 


(P 4 2 

-A—^ r xF 2(q) {x) 


(4.4) 


where | P) now has a mass M and m is the bare qnark mass. 

In the case of quark contributions, the second term in the expression for the bilocal 
current in Eq. (|2.9| ) vanishes. First we evaluate the contribution M 2 given in Eq. (44). We 
obtain, 


M 2 = 


(P^ 

: Q 2 


9- 


5(1 - x) + 


d 2 k 


3 m 2 


¥£kj + (l-x) 

[m 2 (l — x ) 2 + k']_] 2 


r o + (1 — y) 3 m 2 , 

S(l-x)jdydk ± p (1 _ w - fci|2 ) 


(4.5) 


where C f = for SU(N). 

Here we have presented the result without working out the transverse integration to 
maintain a greater degree of transparency. 

The contribution from Mi is split into four parts with additional contributions coming 
from quark mass terms and can be written as follows. 


Mi— j dy e i p+y X (P | i> + \y )( - (d- 1 ) 2 + m 2 )^ + (0) | P) 

+ g—pr, 7 [ dy~e~^ p+v x (P \ i/i + \y~) (id ± .a ± + 7 °m)cr L . 7 b L ( 0 )';/’" t ~( 0 ) | P) 

7TQ- J 

+ g— [dy~e~ 2 P+v X (P I i/j + \y~)a ± .A ± (y)(id ± .a ± + I P) 

nQ z J 

+ g2 ~ql J dy~e~^ r+y ~ x (P I i/j + \y-)A ± (y).A ± (0)'i/; + (0) \ P) (4.6) 

= M a i + M b , + M\ + Mi (4.7) 
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Since the operators in Eq. ( |4.3|) are taken to be normal ordered, the contribution of 
vanishes to order g 2 . 

Explicit calculation leads to the diagonal Fock basis contributions 


(M 


1 Jdiag 


= Mi = 4 


(P±) 2 

Q 2 


-x~ 


4 < 1 - x)+ h°f 


Y^-k']_ + (1 — x) 3 m 2 
[m 2 ( 1 — x ) 2 + k 2 ] 2 

T^ k ± + (! - yfm 2 , 

[m 2 (l — y) 2 + A: 2 ] 2 ' 

l T^r k ± + (1 - y) 3 m 2 


d 2 k i 


— 5(1 — x) I dyd 2 k± 




1 -c, 


8^ r 3 


dyd 2 k 


[m 2 ( 1 — y) 2 + A: 2 ] 2 


4 Cil_ 

Q 2 87t 3 


fci + (i - .t) 


d 2 k±(k~ L + m 2 ) 


3 m 2 


[m 2 ( 1 — x) 2 + A; 2 p 


(4.8) 


The first term here explicitly cancels the term given in Eq. 
Off-diagonal contributions 

^ + 


- / d 2 k j 


A; 2 + m 2 (l — x) s 


(1 — x)[m 2 (l — x) 2 + Af[]- 


Adding all the contributions, we have, 

4 C f g 2 


F t lT,)(x) = ^(1 -x)+ ^’ 


Q 2 


Q 2 87T 3 L 
2 


d k±(k ± + nr) 


Q 2 7T 3 L 


^fc 2 +(l-x) 3 m 2 „ f , ^fA ; 2 +(l-y) 3 m 2 

x ku-*)» + 4]» - ^ ^ 

, 2 A : 2 + m 2 (l — x ) 2 


[m 2 (l — y) 2 + A; 2 ] 2 - 


dk 


— 5(1 — x) / ch/d 2 A;_L- 


m 2 (l — y) 


(4.9) 


(4.10) 


(1 — x)[m 2 (l — x) 2 + A: 2 ] J [m 2 (l — y) 2 + A; 2 ]-" 

Here we have used M = m, since the difference that it entails is higher order in the coupling. 
Note that we are getting back the free quark answer once we switch off the interaction. Also, 
the dressed mass-less quark answer can be easily regenerated by putting M — m — 0. Note 
that the A;x-integration now produces logarithmic divergences with the expected quadratic 
ones, as we remarked earlier. 

To check the sum rule explicitly, we evaluate the RHS of Eq. ( |3.3| ) next. A straightfor¬ 
ward evaluation leads to 
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(P I 0j~(o) I P) nondiag = -C f i / dxd 2 k k2± + m2(1 x)3 


2-7T 3 


x(l — x) [m 2 (l — x ) 2 + A 2 ] 


(4.11) 


(P I C"(°) I P)*w - I », ++ ( 0 ) I P> 4 «» = 


2m 2 + 2 [ dxd 2 k± 
87 r 3 J 


(P 

k? ± + (l-x)m 2 ^k? ± + (l-x) 


3 m 2 


x 


(4.12) 


[m 2 (l — x) 2 + A: 2 ] 2 

Adding the diagonal and off-diagonal contributions from the fermionic part of the Hamilto¬ 
nian density and multiplying it by ^ one obtains the RHS of the sum rule. Comparing it 
with the integral of ——, where Fl is given in Eq. ([4.10D , one easily sees that the sum rule 
is verified. 

To see the connection of F L with the fermionic mass shift, we calculate the contribution 
of the gluonic part of the energy momentum tensor 9 + ~ to the sum rule for the total Fl. 
Explicit calculation gives, 

(1 + x)k\ 1 


(P | ( 0 ) | P) n0 ndiag = ~ C f^ J dxd 2 k A 


(1 — x) 2 [m 2 (l — x) 2 + A 2 ] 


(4.13) 


( p ±)2 

(p 1 e+-(o) 1 p) iim - AA<p 1 c(°) 1 pw» = 


{p*f 


q 2 r k 2 , AT2_A^ + (1 — x) 3 m 2 

20f 8 ^ / dxd k± (1-x) [m 2 (l - x) 2 + A 2 ] 2 ' (4 ' 14) 


Thus, we get, 


2 , (l±al k i + (i, :X ) 3 

I dxd 2 k \ ———— X v ' 


(■P I ^ +_ (°) + 0 J"(°) I P)nondiag = 

1 


m 


2tt 3 


x(l — x) 


[m 2 (l — x) 2 + A 2 ] 


(4.15) 


(P ± ) : 


(P I <>*-(0) + «p( 0 ) | P) dia! - TpLy{P I « ++ + « 9 ++ ( 0 ) I P)*« = 

1 


o 2 r , + (1 - n 3 m. : 

2C/44 / (toft 1_r 1 


87 r 3 J x(l — x) [m 2 (l — x) 2 + A 2 ] 

Adding diagonal and off-diagonal contributions, we get, 


(4.16) 
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(4.17) 


(P | 0 + '(O) | P) - | 0 ++ (O) | P) = 


c 


/: 


9 


dxcPk 


(p + y 

+ (1 - xfm 2 


47r 3 J " L x(l — x) [?tx 2 (1 — x) 2 + A: 2 ] 

Note that this result is connected to the full fermion mass shift 5p{ in second order 


perturbation theory. We have (see Eq. (4.10)) in Ref. |T 8 

{ -^kl + (1 - x) 


Sp i = 


J_ r A _ 

2 P+ / 4tt 3 


dxd k I 


3 m 2 


x(l — x) 


[m 2 (l — x) 2 + k jJ 


(4.18) 


B. Dressed Gluon 


In this section we check the sum rule explicitly for a dressed gluon target. We consider 
the gluon to be composed of a bare gluon and a quark anti-quark pair. 


i»=<M f (P, A) I 0) 


E 




dk2 ^ ^2(2it) 3 P+5 3 (P - h - k 2 ) 


<71,02 


x /2(2tt) 3 A:+ J y / 2(EjEf 
02 (P,C 7 I A:i,cti; fc 2 ,o- 2 ) 6 t (/,’i,ai)d t (/c 2 ,cr 2 ) | 0). 


(4.19) 


The target gluon and the bare quark and anti-quark masses are taken to be zero. Note 
that, to the order g 2 , there will be a contribution from the two-gluon Fock sector due to the 
non-abelian nature of the gauge coupling. For simplicity, we exclude that contribution. It is 
easy to incorporate that contribution by trivially extending our calculation presented here. 

Fl can be written in terms of A4 1 and given in Eqs. ( |4.3| - |4.4|) , where | P) now stands 
for the dressed gluon represented by Eq. (|4.19|). Explicit calculation gives, 


>±112 


Mo = 


a (P ) j-idressed—qluon 

4 ~ 0 ^ xF m 


X 2 (pJ-)2 g2 

Q 2 ^ 


NfTf[x 2 + (1 — x) 2 ]/nA 2 


(4.20) 


Here Tf = | and Nf is the number of flavors. 

Mi is again divided into four parts as in Eq. ( jf.7| ) and explicit calculation in this case 
gives the following. 
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(4.21) 


Adi (diag) Adi (a) 


x 2 (P L ) 2 g 2 

g 2 ^ 

A 2 S 2 


A^T^x 2 + (1 — x) 2 ]/nA 2 


+ + (1 - x ) ] 


g 2 7r 2 


Adi (off—diag) Ad 1 (b) A Ad 1 (c) 

- _ 

g 2 7r 


2_2^/-^(l X ) 


(4.22) 


Thus, we get, 


A 2 


U = ^JV/T,+ (1 - x f - 2(1 - x)) 


Q 


n 


(4.23) 


On the other hand, we get, 

(PU 2 


< p i «r<°) i p)«<* - i ^ + (°> i p w= A 2 N f T 'hi J 


and 


(P | », + -(0) | P)o„-^ s = -A 2 N,T, 


X 2 + (1 — x) 5 

x(l — x) 


(4.24) 


0 


dx[ 


x 2 + (1 — x ) 2 


2 rc 2 J x(l — x) 

Adding diagonal and off-diagonal contributions, we get, 

(P 2- ) 2 / 7 n I /)++/n\ I n\ A 2 at m ^ f j + (1 “ x) 2 


(4.25) 


(P | 0+-(O) | P ) - }^<P I ^ ++ (°) I P ) = -A 2 ^/^ 


(4.26) 


(P +) 2V 19 v /l 7 p p 4yr 2 y 1 x(l-x) r 

Note that this result is connected to the gluonic mass shift 5q 2 due to pair production, 
since the contribution from the gluonic part of the energy-momentum tensor 6 ^~ in this 
case vanishes. In the massless limit, we have (see Eq. (4.40) in Ref. ||18||), 

1 . o „ „ o 2 f . r x 2 + (1 — x ) 2 


5q 2 = 


-P NtT 'h / dx[: 


(4.27) 


2P+ J J 47t 2 J " L x(l — x) 

4 From Eq. ( 4.23 ) we compute / dx^. Since x-integration is from 0 to 1, it can be written 
in the following form. 


, f l A g 

dx ~ = ~W‘ '' w 


2 “ r x 2 + (1 — x ) 2 

dx[- 1 


(4.28) 


Comparing Eq. ([i.26| ) and Eq. 
target. 


x(l — x) 

, one explicitly verifies the sum rule for a dressed gluon 
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As we have emphasized, in the bare theory, the twist four longitudinal structure function 
is afflicted with divergences. We have to add counterterms to carry out the renormalization 
procedure so that we have physical answers. The sum rule for the bare theory clearly shows 
that the quadratic divergences generated are directly related to the gluon mass shift in 
second order light-front perturbation theory arising from an intermediate quark anti-quark 
pair. In order to ensure a massless gluon in second order perturbation theory, we have to add 
the negative of the shift as a counterterm. After adding the counterterm, the gluon mass 
shift in second order perturbation theory is zero and the twist four longitudinal structure 
function for a massless gluon becomes zero. Thus, after renormalization, the sum rule is 
satisfied, with a trivial (i.e., zero) gluon longitudinal structure function. 


V. 1+1 DIMENSIONAL QCD: EXPLICIT CALCULATIONS 


In this section, we turn to two-dimensional QCD in order to test the sum rule given in 
Eq. (|3.4J) explicitly in a non-perturbative context. In 1+1 dimensions, in A + = 0 gauge, we 
have, 


rl dr 2 r ~\ 

/„ - F m( x ) = Q2< p [»r(0) + + (0)] I A 


(5.1) 


with 6 +~ = 2m 2, 0 + ^^+'0 + and = —4g 2 7/i + ^r a, 0 + ^ +t T a, 0 + . We consider the stan¬ 

dard one pair (qq) approximation to the meson ground state. Explicit evaluations show 
that 


F ko) 4 ,* f \ , ( \ 

x Q z x(l — x) 


(5.2) 


and 


f d /J5L = A_ )c / r d q iynx) m^m 

Jo x Q 7r jo jo [x — y) z 

where 'ijj(x) is the ground state wave function for the meson. Thus 


(5.3) 


rl dx 


I o x 


nr=4/ 


el 


L ( X ) = q5J 0 dxip*(x) 


rn 


9 


y) - 


I n M x ) ~ Cf—^*{x) / dy- 
L o:(l — x) 7 r jo (+ — y) z 


(5.4) 
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By virtue of the bound state equation (’t Hooft equation) obeyed by the ground state wave 
function iJj(x) for the meson 


= ~ Cf ^ J dy 


V>( y ) - ^{x) 

(x-y) 2 


(5.5) 


together with the normalization condition Jq 1 dxip 2 (x) = 1, we easily verify that the twist 
four longitudinal structure function of the meson obeys the sum rule 


rl dx 


r>T=4 


= ™(-P I « +_ (0) I P) = 4 


M 2 


(5.6) 


Jo x - Q 2 ' ' ' Q 2 

In the same model, the contribution to the twist two structure function from the fermionic 
constituents is given by 


F 2 ( q )(x) — (x + 1 — x) / tp*(x)^(x) = ip*(x)ip(x). (5.7) 

Note that, since there are no partonic gluons or sea in this model, the longitudinal momentum 
of the meson is carried entirely by the valence quark and anti-quark. Thus the momentum 
sum rule is saturated by the fermionic part of the longitudinal momentum density. On the 
other hand, light-front energy density is shared between fermionic and gauge bosonic parts 
and as a consequence the fermions carry only a fraction of the hadron mass. This seemingly 
paradoxical situation further illuminates the difference between the physical content of the 
F 2 and Fl =A sum rules. 

To get a quantitative picture, next, we explicitly calculate the structure functions F 2 ( q )(x) 
and for the ground state meson in two-dimensional QCD. We have parameterized 

the ground state wave function as ip(x) = Afx s (l — x) s and determined the value of s 
variationally by minimizing M 2 for given values of m 2 and g 2 . The factor M is determined 
from the normalization condition Jq 1 dxip*(x)'ijj(x) = 1. The resulting structure functions are 
presented in Fig. 1 for two different values of m 2 . 

Since both the quark and anti-quark have equal mass in the model, both structure 
functions are symmetric about x — When the fermions are heavy (Fig. 1(a)), the system is 
essentially non-relativistic and the structure functions are significant only near the region x = 
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When the fermions become lighter (Fig 1(b)), contribution to the structure function from 
the end-point regions become significant indicating substantial high momentum components 


in the ground state wave function. Note that —^ measures the fermion kinetic energy (in 
light-front coordinates). The exponent s in the wave function is a function of the fermion 
mass and s decreases as m decreases. In the massless limit, s vanishes |19|] so that the 
wave function for the ground state becomes 'ip(x) = 9(x)9( 1 — x). This results in a flat 
F 2 structure function. However, because of the presence of m 2 , FW? vanishes. Because of 
an exact cancellation between the self-energy and gluon exchange contributions, the gluonic 
part of the FJ =4 also vanishes. Thus the sum rule is satisfied exactly since, in the zero quark 
mass limit, the ground state meson is massless in two-dimensional QCD. 


VI. PARTITION OF THE HADRON MASS IN QCD 


As is well-known, experiments that measure the twist-two part of the F 2 structure func¬ 
tion yield information on the fraction of longitudinal momenta carried by the charged parton 
constituents of the hadron (quarks and anti-quarks). The sum rule we have derived yields 
other useful information about the hadron structure. Namely, our sum rule shows that ex¬ 
periments to measure the twist four part of the longitudinal structure function will directly 
reveal the fraction of the hadron mass carried by charged parton components of the hadron. 
The light-front Hamiltonian provides theoretical insight into this fraction as follows. 

According to our analysis, the twist four part of the longitudinal structure function is 
directly related to the fermionic part of the light-front QCD Hamiltonian density 9 4 in the 
gauge A+ = 0. Explicitly we have 


9+ = a ± .(id L + gA L ) + 7 °m -— a ± .(id ± + gA L ) + 7 °m ^ 


id + L 


( 6 . 1 ) 


Thus we have the fermion kinetic energy contribution given by 


7<7r,.) = 2 F\ ~ (9 1 ) 2 + ro 2 ]V> + 


( 6 . 2 ) 


and the interaction dependent part given by 
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(6.3) 


d t(int)= 


+t 


a L .A^ .id L + 7 °m) + (a ± .id ± + '-y°7n)^—a ± .A 1 - ^ 

io + io + J 

A ± 


+ 2 g 2 7/j + 'a ± .A ± a ± .A ± ^ 


Note that the fermion kinetic energy constitutes only a part of the total contribution from 
fermions. Any theoretical estimate of the fermionic part of the longitudinal structure func¬ 
tion necessarily has to involve off-diagonal contributions from Fock states differing in the 
number of gluons by one and two. 

It is important to emphasize the difference between equal time and light-front Hamil¬ 
tonians in the context of our calculations. The equal-time Hamiltonian contains the scalar 
density term accompanying the quark mass m. I 11 contrast, the quark mass appears 
quadratically in the free part of the light-front Hamiltonian. Recently the question of the 
partition of hadron masses in QCD has been addressed by Ji [^tj in the context of the equal¬ 
time Hamiltonian and in terms of twist-two and twist-three observables. In his analysis, the 
extraction of the fraction of the hadron mass carried by the fermion constituents is not 
straightforward because of the presence of the scalar density term. The hadron expectation 
value of the strange quark scalar density remains unknown (experimentally). Our analysis, 
however, shows that the twist four longitudinal structure function, once extracted experi¬ 
mentally, directly yields the fraction of the hadron mass carried by fermionic constituents. 


VII. DISCUSSION AND CONCLUSIONS 

To gain physical intuition on the twist four longitudinal structure function and to un¬ 
derstand the occurrence of quadratic divergences and the associated renormalization issues, 
we have studied the twist four longitudinal structure function in an approach based on Fock 
space expansion methods in light-front field theory. We have identified the integral of ^ 
times the twist four part of the fermionic contribution to longitudinal structure function as 
the hadron matrix element of the fermionic part of the light-front QCD Hamiltonian density 
in the light-front gauge apart from an overall constant. We have tested this relation to order 
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g 2 in QCD perturbation theory for both dressed quark and gluon targets. Our result shows 


that quadratic and logarithmic divergences in the twist four longitudinal structure function 
are directly related to mass corrections in the light-front theory. 

By investigating the mixing of operators in the flavor singlet channel, we have recently 
derived [|T^] a new sum rule which involves the invariant mass of the hadron. The validity of 
the sum rule has been explicitly checked in two dimensional QCD (’t Hooft model). To get a 
qualitative picture of the twist four structure function we have computed numerically both 
F 2 and Fl structure functions in the ’t Hooft model using the ground state wave function 
calculated using a variational ansatz. 

We have also discussed the implication of our results for the problem of breakup of 
hadron masses in QCD in terms of fermionic and bosonic constituents. We have emphasized 
the differences between equal-time and light-front formulations relevant for this study. 

Our results indicate that the experiments to measure the twist four longitudinal structure 
function reveal the fraction of the hadron mass carried by the charged parton components. 
Thus these experiments play a complementary role to the longitudinal momentum and he- 
licity distribution information obtained at the twist two level. It is of interest to investigate 
the feasibility of the direct measurement of the twist four gluon structure function in high 
energy experiments. Recent work of Qiu, Sterman and collaborators have shown that semi- 
inclusive single jet production in deep inelastic scattering |H[] and direct photon production 
in hadron nucleus scattering [^] provide direct measurement of twist four gluon matrix 
elements. 

On the theoretical side we note that some significant progress has been made recently in 
the bound state problem in light-front QCD |23j based on similarity renormalization group 
method. In the near future, we plan to undertake a non-perturbative calculation (utiliz¬ 
ing Fock space expansion and Hamiltonian renormalization techniques) of the longitudinal 
structure function for a meson-like bound state. Such a calculation will undoubtedly help 
to check the validity of current phenomenological models [B| based on simple assumptions 
2 4j employed in analyzing the data. 


20 











Another important problem is the scale evolution of the twist four structure function 
which is far more complicated than the twist two case. Recently we have provided a physical 
picture of scale evolution of the F 2 structure function of a composite system in terms of multi- 
parton wave functions in momentum space [0. We plan to carry out a similar analysis of 
the twist four longitudinal structure function elucidating all possible scale dependencies and 
their physical interpretation. 
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APPENDIX A: SUMMARY OF NOTATIONS AND CONVENTIONS 


The hadron tensor relevant to unpolarized electron-hadron deep inelastic scattering is 
given by 

= ( - f r + )Wi(x, Q 2 ) + (p 11 - [P v - £jL<?)w 2 {x, Q 2 ). (Al) 

The dimensionless functions 


and 


Fl(x,Q 2 ) 


2 


Wi + [M 2 




(A2) 


F 2 (x,Q 2 ) = vW 2 {x,Q 2 ) (A3) 

are the unpolarized structure functions. 

We have defined, — Q 2 = q 2 = q + q~ — (q L ) 2 ■ 

The light-front coordinates are defined by x ± = x° ± x 3 . 
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The constraint equation for the fermion field which follows from the Dirac equation, in 
A + = 0 gauge is given by 

^~( z ) = -^f dy~e(z~ - 2r)[or L .(id i - + gA^) + 7 0 m]^ + (y“), (A 4 ) 

where the antisymmetric step function 

e(x~) = --P f—e^ x ~. (A 5 ) 

7 T J CU 
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j?t=4 

1. Fermionic contributions to the structure functions F 2 (x) and —— for the ground state 

meson in the’t Hooft model for two different values of m, the quark mass, (a) m = 5, 

s = 4.96. (b) m — 1, s = .70. The parameter s appearing in the wave function is 

2 

determined by a variational calculation. We have set Cf — = 1. 
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